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SUMMARY 


This note presents the time-dependent, turbulent mean-flow, Reynolds 
stress, and heat flux equations in mass-averaged dependent variables. These 
equations are given in conservative form for both generalized orthogonal and 
axisymmetric coordinates. For the case of small viscosity and thermal con- 
ductivity fluctuations, these equations are considerably simpler than the 
general Reynolds system of dependent variables for .. compressible fluid and 
permit a more direct extension of low speed turbulence modeling to computer 
codes describing high speed turbulence fields. 



INTRODUCTION 


In recent years, significant advancements in both computation techniques 
and computer hardware have led to sophisticated prediction techniques for 
describing complex compressible fluid flows. The most general techniques avoid 
boundary- layer approximations and solve the complete time-dependent, conservation 
equations of mass, momentum and thermal energy, and possibly species. (See 
ref (1) and (2).) The transport mechanisms that have been employed in these 
techniques are consistent with molecular 01 laminar processes. Methods such 
as those presented by MacCormack or Carter have been shown to apply to a flow 
as complex as a shock-wave boundary- layer interaction, even when the adverse 
pressure gradient generated by the shock wave is strong enough to cause separa- 
tion of the boundary- layer flow. 

The successful methods for laminar flows have not, however, been carried 
over to the computation of turbulent flows — particularly strongly interacting 
flows. Turbulent-flow solutions, which are presumably contained in the instan- 
taneous viscous flow conservation equations, cannot be obtained at present 
because the small grid size required for accurately resolving the decay pro- 
cesses of turbulence cannot be accommodated in present computers. Therefore, 
these instantaneous equations must be time averaged to obtain a system of 
equations with spatial characteristics consistent with present computers. In 
the time averaging process there is a loss of phase and wave-length information 
of the turbulent motions. Although this information may be critical to the 
ultimate understanding of turbulence, time averaging has yielded most of the 
significant results related to the mean flow in boundary layers and may prove 
equally valuable for the more complex flow fields of current interest. 



The first step i establishing time -averaged turbulent flow equations is 
the expansion of all the dependent variables into the sum of two terms, one 
slowly varying and the other randomly varying at the high frequencies charac- 
teristic of turbulence. With compressible flows there is some freedom in the 
choice of how this expansion is accomplished. For example, instantaneous 
velocity components can be expanded either as the sum of two velocities, as in 
the classical Reynolds expansion, or as the sum of two mass fluxes divided oy 
the mean density, the mass averaged dependent variable expansion, first suggested 
by Favre (ref. 3). As an example of Reynolds 1 expansion, Donaldson (ref. 4) has 
written the complete system of equations for compressible flow in generalized 
tensor form so that they can be particularized to any desired coordinate system. 
Donaldson's equations explicitly contain density fluctuations as well as velocity 
fluctuations, thus, adding considerable complexity as compared with the constant 
property case. This added complexity can be reduced, as shown in ref. 5 and 6, 
if the expansion of the dependent variables is carried out using mass averaged 
expansions. These equations essentially exhibit a term by term correspondence 
with the constant property case. This correspondence can be important if the 
basic turbulence mechanisms in compressible flows remain similar to those in 
constant property flows (the so-called Morkovin hypothesis), thereby permitting 
convenient extension of the large body of experience existing with incompres- 
sible flows. 

The work of refs. 5 and 6 is limited, however, in that the Reynolds stress 
and turbulent heat flux equations were not presented and, further, the equations 
were confined to the steady state and expressed in cartesian coordinates. The 
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general objective of this note is to remove these limitations by 


1. adding the time-dependent terms to all the equations 

2. presenting the Reynolds stress and heat flux equations 

3. employing a generalized orthogonal coordinate system 

4. casting the equations in conservative form, consistent with existing 
computer codes (refs. 1 and 2) 

5. particularizing the equations to cylindrical coordinates for com- 
parison with axisymmetric flow data (e.g. ref. 7). 
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SYMBOLS 


total internal energy per unit of volume, ph - p + p 
specific enthalpy 

u.u. 

total specific enthalpy, h + 
static pressure 

period of integration for turbulence quantities 

molecular Prandtl number 

heat flux in j th direction, - A |^- or - 

J r dx i 


u- u- 
J J 


dummy variable of integration 
time 

temperature 

velocity component in i th direction 
molecular thermal conductivity 
molecular viscosity 
density 

instantaneous shear stress in i th direction on a surface normal to the 



2 . V 

3 ij 9x k 


j th direction. 



TURBULENCE EQUATIONS IN ORTHOGONAL TENSOR FORM 
In the derivation of the mean flow and turbulence equations, orthogonal 
tensor notation will be employed. The usuil conventions of a repeated sub- 
script indicating summation over the entire range of the indices and a comma 
representing partial differentiation will be used to make the equations compact. 
Orthogonal tensors were adopted to permit working directly with the physical 
components, thereby avoiding the need for identifying specific tensors as 
being covariant or contravariant. The use of orthogonal coordinate systems 
does not cause significant loss of generality in most fluid mechanical 
applications . 

The time-dependent, conservation equations for a compressible fluid with 
zero external body forces are 
Continuity 

P t + (puj . = 0 (1) 

* u J y J 


Momentum in the ith direction (i = 1, 2, or 3) 


(pu.) _ + (pu.u. + 6..p - t..) . ~ 0 
v i',t i 3 l 3 ij >3 


( 2 ) 


Thermal Energy 


( p H * P\ t + (P u j H + " Vij\j = 0 


(3) 


to convert tnese equations to time averaged equations, the following expansions 
of the dependent variables are employed: 

p * p p ' 


T. . 

13 


T. . + T. . 

13 13 


q. + q. 
j S 3 


( 4 ) 
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and the mass averaged dependent variables 

n 

u. = a. + u. 

1X1 

H = H + H" (5) 

h = h + h" 

The superscripts - , - represent time averaging over a period P that is long 
compared to the fluctuation period of the turbulence but short compared to 
the time constant for changes within the flow field as a whole. The super- 
script ' or " represents the instantaneous fluctuation about the corresponding 
time averaged quantity. Thus, we set 

J 0t +(P/2) 

I ° A(t)dt (6) 

t 0 -(P/2) 

then. 


3A(t o ) 

at 




V (P/2) SACi dt 

. at 

(P/2) 


ft A <V 


For example, from Eqs. (4) and (6) 


1 |V (P/2) . 

p(t ) - p I (p(t ) . p (t) jclt 

0 P Jt o -(P/2) 


(7) 


( 8 ) 


a P(t Q ) + P'(t 0 ) 


(9) 


From Eq. (9), 


p' * C 


( 10 ) 
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Thus, in general 


• ■ H: 


t Q + (P/2) 


A'(t)dt = 0 


( 11 ) 


(P/2) 


a single fluctuating quantity raised to the first power averages to zero. 

The mass-averaged dependent variables are defined (see Appendix B) such 
that the mean velocity is 


0 . 


1 



and the fluctuating velocity averages to 



From Eqs. (4), (5) and (13) 



so that uV is not zero in contrast to Eq. (11). 

Similarly, the mass averaged static enthalpy is 


( 12 ) 


(13) 


(14) 


and 



ph" = 0 


Again h" is not zero, but rather 


h" = - 


P 


( 15 ) 


( 16 ) 
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Mean Flow Equations 


When the instantaneous density and velocity, 
are inserted into the continuity equation (1) and 
time averaged according to Eq. (6), there results 

p t + (pa.) • = 0 


express as Eqs. '4) and (5), 
the resulting equation is 


(17) 


This time-averaged continuity equation has identical form to Eq. (1) except 
that the instantaneous quantities p and u are replaced by the time averaged 
quantities p an> u. 

Similar operations of expanding the dependent variables according to 
Eqs. (4) and (5) and then Reynolds averaging converts the momentum equations, 
Eq. (2), and the energy equation, Eq. (3), to 


and 


Cpa i>.t * 


[ 5 Vi 


<5. .p 


- (t . . 

\ 13 


pu'.’u 

1 




(18) 


(pH - p) 




pu.H 

J 


pu'.'h" - u. ( i. . 
3 i \ 13 


/ pu'.'u'.' \ 

- pu'.'u'.' ) - u'.' (t. . } 

13/ i \ 13 2 /J ,3 


(19) 


Again, their form is generally identical with their laminar flow counterparts 
except for the addition of the Reynolds stress, heat flux, and the additional 
mean energy dissipation terms given below: 


Reynolds stress 



pu'.'u'.' 
i 3 


( 20 ) 
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( 21 ) 


Reynolds heat flux = q^ = puVh" 


/ PuVu'W 

Mean energy dissipation = uV It.. 2~j 


Turbulence Growth Equations 

The growth of the Reynolds stresses are found by manipulation of Eqs. (1) 
and (2) into a moment of momentum equation 


(pu iV,t * (eu iVj>,j * v.i * v,k • vijj - u i T kJ ,j ■ 0 


This equation is then expanded according to Eqs. (4) and (5) and time averaged 
according to Eq. (6) to yield the sum of the mean flow terms 

KV.t * hW.i 

■ -v.i - V.k * “k ( ; ij - ™H),j * “i( ; kj - ^).j (24) 

and the turbulence quantities (growth of Reynolds stresses) 


( pu i u k),t * (v u i u k),j 


- ("tl) 4 !.] - ( PU i U k U j ) ,j 

- (v).i - (v).k * p (It * U i.k) 

* (V7j),j * (vTj).j - V '.j • VL 





Eq. (24) follows directly from Eqs. (17) and (18); so that the separation of 
Eq. (24) and Eq. (25) as indicated is not arbitrary. Eq. (25) corresponds term 
by term to the constant property equation in Ref. 8, except with p occurring 
inside of the operator signs of time-averaging and differentiation. 

The turbulence kinetic energy equation, with 


k = 


k pu i u i 


(26) 


and 


k 



(27) 


is 


(pk) + (pu.k) 

* J > J 


* -( pu i u j)“i,j - ( pu j k )j - ( u i p ),i * pu i,i * ( u i T ij)j 


T. .UV 
1 ) 


(28) 


The terms on the right are, sequentially, the turbulence production, the 
turbulent diffusion, the pressure diffusion, pressure work, molecular diffu- 
sion, and dissipation of the turbulence kinetic energy. Eq. (28) is consistent 
with Rotta's and Townsend's (Ref. 9) constant property equations. 

The equation for the growth of Reynolds heat flux, Eq. (21), follows from 
the energy equation, Eq. (3), combined with Eq. (23) with i - k and the 
equation resulting from h times Eq. (?) 


< pu i h) , t * (pu iV : ,j ■ v.t - w,k * h >\i 


+ u.q. . - hr.. . - u . t. . u. . = 0 
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(29) 



Again, when this equation is expanded according to Eqs. (4) and (S) , and the 
resulting equation is time averaged according to Eq. (6), there is obtained the 
sum of the mean quantities 


■ v.t * "i * pu j h " - Vkj - u k T kj],j 

- “iV.k - a i“i? 7 k * “iVkj.i * WkJ.j * fi P,i - - ™H).j * 0 


(30) 


and Reynolds heat flux terms 


( ou P"),t * ( a j pu i h ").j * ( ou j h ”) a i, j * ( pu i' u i') h .j * ( pu j u i h "),j - “P.t 


V?,k - “i-'iP.k * h "P,i - h " T ij,j * “Pj.j - V; u i T ki ' “kjVkj 


(31) 


The proper separation of Eqs. (30) and (31) follows from the use of moments of 
Eqs. (18) and (19) to yield Eq. (30). 

THE AXISYM4ETRIC-FLOW EQUATIONS 

The equations given above in orthogonal tensor form are particularized in 
this section to an axisymmetric system (i.e. a cylindrical coordinate system 
with * 0 and w = 0). These equations are obtained using the particularized 

O v 

expansions and coordinate nomenclature given in the Appendix. 

Mean Flow Equations 

The mean flow equations Eq. (17) through (19) particularize to 
Continuity : 

illEl + + Hrpv) = o (32) 

3t 3x 3r 
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x-momentum 


^-(rpQ) + ^-[r(paQ + p - r a ♦ pu"u")] ♦ ^r[r(p\Kl - T xr + pu"v")] = 0 (33) 

r-momentum 

-^(rptf) + ^-[rCpa^ - x^ x ♦ pu"v")] ♦ ^[r(p« + P - T rr + pv"v")J 


+ Tqq - p - pw"w" = 0 (34) 


Total thermal energy 


* h r [‘ 


a(ft + p) + q x ♦ pu"h" - Q(t xx - pu"u") - V(T rx - pv”u”) - u"(t xx - ) 

" ( T rx - e}f ¥ L ) - w "( T ex - £S r-)] * ^{*< s * rt * s * 


Q(x xr - pu'V') - 9(1^ - pv’ 


•V) - u"( V - *F) - v( v - <*?) 



(35) 


These equations may be written in a condensed vector form convenient for numerical 
computations (ref. 1 or 2) as 


9rU 3rf 3r$ ^ | 

ir + ir + ir** 


0 


(36) 


whers, for example. 



( 37 ) 
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These equations may also be written in an alternative form by dividing through 
by r to give 


3U 

3t 


+ 



3r r v 


+ A) = 0 


(38) 


We may combine the two inhomogeneous terms as 

I .«*! 

r 


In this form ^ becomes the following vector: 



In the alternative form of the equations, the transformation from the axisym- 
metric case to the planar flow case is simple in that the vector B vanishes 
and the remainder of the equation is unchanged, i.e., we have for the planar 
case 


jW + + il - 

3t 3x 3r “ 


(40) 
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Turbulence Growth Equations 

To affect a closure of the above equations (in either of the two forms 
presented) we require equations for pu"^, pu"v", pv" 2 , pw" 2 , pu"h", pv"h" as 
well as the terms in the thermal energy equation which have the form: 



i puVuv) 


Relative to these latter terms, Laufer (ref. 5) points out that, although all 
previous analyses have neglected them, they can be of the same order as the 
turbulence production terms . 
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Reynolds stress equations in the axisymmetric form are 



where i and k take on values from 1 to 



The equation for the turbulence kinetic energy in axi symmetric form is 


& J .. 




equations for the turbulent heat flux in axisymmetric form are 


i 

I 


I 

\ 

I 


£ 


ol u 
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where i takes on values 1 and 2, for heat flux in axial and radial directions, respectively. 



COMMENTS ON THE USE OF MASS AVERAGED VARIABLES 

The use of the mass -averaged variables results in mean flow turbulence 

equations that contain a minimum number of terms, since terms involving p' do 

not appear explicitly, and posses:, almost a term by term correspondence with 

their incompressible counterparts. On the other hand, the complete expressions 

for the molecular shear terms t. . contain extra terms. It will be shown that 

13 

this latter condition does not seriously complicate otters. 

Instantaneously, 


t. . = y mu. . + u. \ 6. .u, ,1 ( 4 ‘ 

13 l\ 1,3 3,1/ 3 13 k,kj 

when the Stokes criterion for the second viscosity coefficient is invoked. For 
compactness, we write the quantity i.i the braces as so that 


x 


ij 



( 45 ) 


When Eq. ( 45 ) is expanded in the Favie variables, there results 


1.1 


= X . . 
11 


T ! . = (y 

13 


y") (S. . + S'.'.) 

J 13 ir 


( 46 ) 


Here the viscosity is identified with the mass averaged enthalpy of Eq, ( 15 ) or 
the equivalent mass averaged temperature. (See Appendix B for the interrela- 
tionship of the Favre viriable with the usual Reynolds expanded quantities.) 
From Eq. ( 46 ) 


x.. * x.. ♦ x.. = yS.. 
ij U 11 il 


. “c" . "2 . Il 0 ” 

+ yS.. + y S.. + y S.. 
13 13 13 


( 47 ) 
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Eq. (47), time averaged, yields 


~x ~ c " , ii; , tt c it 

t.. = yS. . + |iS.. + y S.. + y S.. 
lj ij lj lj ij 


and 


* c t " "n . '( c" c" ^ ^ "c" "o" 

t. . = S..(y - y ) + y(S. . - S. .) + y S. . - y S. . 

lj lj" i] ij" ij ij 


(48) 


(49) 


In customary Reynolds expansions, not the mass-averaged coordinates, time 
averages of first order turbulent quantities vanish, which would have elim- 
inated half the terms in Eqs. (48) and (49). Thus the mass-averaged notation 
appears to complic^L; the expressions involving the molecular viscosity. This 
additional complexity is mitigated through arguments of the relative magnitude 
of the individual terms of Eqs. (48) and (49) when equations as complex as 
Eqs. (41) through (43) are considered for solution. 
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APPENDIX A 

Orthogonal Tensor Relationships 

The following definitions apply for an orthogonal curvalinear coordinate 
system (Ref. 10). The metric, or length of a line element is 


ds 2 = h 2 d<| + h 2 dx 2 + h 2 dx 2 (A-l) 

where hy, h?, and hj are the metric elements of the orthogonal coordinates 
xj , X 2 and X 3 .* 

The g radient of a scalar , a vector quantity, is given by 


i _ 1 

^,i h. 3x. 
x 1 

The gradient of a vector , a tensor quantity, is 

3A. 3h. 

k 


(A-2) 


U h. 3x. + d \ 


A. 3h. 


h.h. 8x. 
1 3 1 


(A- 3) 


where 


6^ (Kronecker 6) = 1 i = j 

= 0 i f j 

Eq. (A-3) can be applied to define the directional derivative of a vector , 
namely 



B.A. . 
3 1.3 


3 


y^ jlIa + y^ 

Z-> h 8*4 Lj 

3 J j 


B.A. 9h. 
i ] 1 

h.h. 3x . 
1 3 3 



B.A. ah. 
_i_J_ — I 
h.h. 9x. 


(A-4) 


♦Note that the summation convention of a repeated index is not used in this 
appendix. 
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This disagrees with Ref. (10) in that the subscripts of A and B in the second 
term on the right are interchanged there. No significant consequence occurs, 
howeve?. , since in Ref. (10) both B and A ultimately represent the same quantity. 


Again, Eq. (A-3) can be summed to yield the divergen ce of a vector A^ 


^ , A j,j hih 2 h 3 [ 


- 4 (hlh,A2) * 15I th > h ^ ) 


] 


The divergence of a tensor can be expressed from Eqs. (A-4) and (A-5) as 


* h71^[^ Ch 2 h 3 T ii) + ^Chihaxiz) + ^_(h 1 h 2 T i3 )J 


1 3h i 

+ x. . r— r— tt— + t 


. 3h. 

1 l 


li h.h 1 3Xj T l 2i h.h 2 ax 2 + T 3 i h.h 3 9 x 3 

. 3h 


1 3h i 


. 3h 
1 3 


1 3h l _ _J _ 

T ll h^ ' T 22 h 2 h i 3x i " T 33 h^ 3x i 


Also required in the orthogonal tensor form of the equations is a 
tensor directional derivative analogous to the directional derivative of 
a vector which can be written as 


ZJ V i’i 2 -j [ h j 9x 3 h i h j 3x j h i h j 3x ij 


(A-5) 


(A-6) 


(A-7) 
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In addition, the divergence of a third order tensor is required, ,'nis form 
may be written as 



V—nTi 3 E... E... 3h. E... 9h. E... ah. E... 3h. 

l_j [h. axj igr 3XJ h^ir 3^ h.h. 3 Xj hjh. 3x t 


J 

E. . 3h. E. . 3h. E. , 3h. E. .. 8h. 

ilk j + 12 k j + 13 k 3 _ ijk j 

h.h 3x, h.h, 3x, h h ax h.h. 3x. 

J 1 1 3 2 7 - 3 3 3 3 3 JJ 


(A-8) 


Axisymmetric Case 

When particularizing to the axisymmetric coordinate system of interest 
here, the above general forms are simplified considerably. For the axisym- 
metric case we take a cylindrical coordinate system defined by 


xj s x h x = 1 
X 2 * r h 2 * 1 
X3 = 6 h3 » r 


We may then write the forms discussed previously as follows: 
Gradient of a Scalar: 


♦ 


t 





♦ -I»i 

\3 r 30 


Gradient of a vector: 



3A i 

3Xj + fi j3 5 i3 



r 


6 . 6 . 

i2 3J 



r 


Divergence of a vecto r 




A 

JL 

r 


(A-9) 


(A- 10) 


(A-U) 


(A-12) 
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Directional derivative of a vector 


E B. 9A. 

ri- — - + 6 . 
h j 9x j 1 


B.A. . 


B A 
3 2 


B A 

6. -L-l 


3 r 12 r 


"Tensor directional derivative" 


E B. . 9A. ». 

— • 

hj 3Xj 13 


6 . 


B, A, 
ks 3 


12 r 


Divergence of 2nd order tensor 


B. . 
13,1 


ri2u.i|, . 

l_, h. 3x. r (i 2 
i 


6. B 
13 23 


6. B 1 

i2 33 


Divergence of 3rd order tensor 


E. . A ^ r + - r 6, E. - 6, E ♦ 6 E - 6 E + 5. ,1 

3 3 k »3 Z J h i 3x i r I k3 132 k2 i33 i3 23k i2 33k l2k> 
4 3 J 


) 


The above forms were used, in conjunction with the conditions, 0=0 

IQ 

39 


and d ' sM * 0 deriving the axisyrmcetric equations presented in the body of 


(A-13) 


fA-14) 


(A-15) 


(A-16) 


this 


note 






APPENDIX B 

The Relationship Between The Reynolds And 
The Mass-Averaged Dependent Variables 
In the past, the Reynolds system of dependent variables has been used in 
both analytical investigations and interpretation of experimental data. Since 
the equations discussed in this report are cast in terms of the mass-averaged 
variables introduced by Favre, it is important ir relating past work to this 
report to know how the variables, particularly the fluctuating variables, are 
related in the two systems. 


We begin by summarizing the definition of the dependent variables in the 
two systems. 


Dependent Variable 

Reynolds 

Mass Averaged 

u 

u + u' 

a + u" 

V 

V ♦ v' 

V + v" 

w 

w + w' 

a + w" 

h 

h + h' 

h + h" 

T 

T + T 1 

f + t" 


u + y' 

u + u" 

P 

p + p' 

p + p' 

P 

p + p« 

p + p' 


The tilde indicates a mass-average defined in terms of the conventional 
time average ^s 



(B-l) 
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We also note that the variables p and p are defined identically in terms of 
mean and fluctuating components in bo«‘h systems. 

To reveal the relationship between the systems we examine the character- 
istics of an arbitrary variable a, i.e. 

a = a + a' = a ♦ a" (B-2) 

xrst, we note that in terms of the mass-averaged quantities 

pa = pa ♦ pa" = pa ♦ pa" (B-3) 


and, thus. 


pa" = 0 


Further, we may evaluate a" from the above, since 


0 ■ pa" = pa" ♦ p'a" 


(B-4) 


(B-5) 


and, thus. 


Note that, although a' * 0 in the Reynolds system, a" 4 0 for the mass- 
averaged system. 

In addition, since 


(B-6) 


a * a ♦ a' * a + a" 


(B-7) 


i 
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we can tine average this to show 


J 


a - 5 * a" 


(B-8) 


and then show that 


a" = a ♦ a' - a 


P* a ‘ 


(B-9) 


Multiplying this last equation by p and tine averaging, i.e. 


0 = pa" = pa' - P ? P * p'a' - p'a" 
P 


we see that 


p'a" 


(B-10) 


and, therefore. 



(B-ll) 


With equation (B-li), we may immediately deduce the relationship between 
turbulent transport of momentum and enthalpy expressed in the two systems 
of dependent variables. For momentum transport in the mass -averaged system, 
we require terms of the form 

i* *< 

pa 6 
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which can be expressed in the Reynolds system as follows: 


i»„" " /" i p 'a' \ /„ i p’fi* \ 

pa 6 • p(. - — ) (B - — ) 


We may now examine some specific examples. Consider the turbulent shear 


stress term pu"v": 



pu"v" * 

pu'v' - 

p 'Up 'V' 

CB-13) 


P 

similarly, the 

normal stress term 

pu" 2 is 

just 





/ \ 2 



pu" 2 

= pu' 2 

Vp'uV 

p 

(B-14) 

The turbulent 

heat transfer terms 

are as 

follows 

F^pMT’ 



pu"h" ■ 

pu'h' - 

(B-15) 


P 

and 






pv"h" - pv'h' - £- ' - v 4 p - ' — (B-16) 
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